Frustrated collisions and unconventional pairing on a quantum superlattice 
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We solve the problem of scattering and binding of two spin- 1/2 fermions on a one-dimensional 
superlattice with a period of twice the lattice spacing analytically. We find the exact bound states 
and the scattering states, consisting of a generalized Bethe ansatz augmented with an extra scat- 
tering product due to "asymptotic" degeneracy. If a Bloch band is doubly occupied, the extra wave 
can be a bound state in the continuum corresponding to a single-particle interband transition. In 
all other cases, it corresponds to a quasi-momentum changing, frustrated collision. 
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a. Introduction. The study of electrons in dis- 
cretized models with a periodicity larger than the lat- 
tice spacing has a broad range of physical applications. 
Among them, a superlattice model for ferroelectric per- 
ovskites has been proposed [li, and its rich phase dia- 
gram qualitatively explained 0. The recent experimen- 
tal advances with ultracold bosonic or fermionic atoms 
loaded in optical lattices provide a neat realization of su- 
perlattice Hamiltonians [3]. These may allow for exper- 
imental observation of fermionic d-wave superfluidity Q 
and magnetic phases of cold spinor gases Q. Moreover, 
one-dimensional electron systems can exhibit dimeriza- 
tion due to Peierls' instability @, thus enlarging the pe- 
riodicity of the system, by a factor of two, exactly 
There is growing interest in the theoretical and numeri- 
cal understanding of the quantum phases of superlattice 
and dimerized models, as attested by the attention re- 
ceived in the recent literature . Although essentially 
exact numerical calculations in one dimension can be per- 
formed, much of the physics is hidden in the numerical 
data, and no unambiguous physical interpretation can be 
given in this way. In particular, the exact origin of degen- 
eracy in dimerized models is not known, and may only 
be determined via an exact analytical solution. 

In this Communication, we consider the problem of 
two-electron scattering and binding on a one-dimensional 
superlattice with a period of twice the lattice spacing. Af- 
ter introducing two different Hubbard models with non- 
trivial periodicity, we show that the Schrodinger equation 
can be solved exactly. To do so, we generalize the two- 
body Bethe ansatz [13] of the Hubbard model [ll| in a 
non-trivial manner, and center our discussion on one of 
the models, for concreteness. We find that, surprisingly, 
the collision of two electrons with opposite spins either is 
frustrated — the scattering states are degenerate — or 
necessarily a bound state in the continuum is produced 
upon collision. We also find all (up to four) bound states 
of the system analytically. Our results constitute the first 
quantum scattering problem on a non-trivial periodic po- 
tential, and with both particles being mobile, to be solved 
analytically, and are also relevant to the study of dilute 
quantum gases in optical superlattices where two-body 
processes are dominant. 



h. Periodic Hubbard models. We start with the 
following general Hubbard Hamiltonian for spin- 1/2 
fermions 
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where 



(1) 

j(T (^jo-) is the creation (annihilation) operator of 

a fermion with spin a =^,-1 at site j, nj„ — c^jcr^j<y is 
the number operator, Jj is the tunneling rate between 
adjacent sites j and j + 1, Xj is a single-particle potential 
and U is the on-site interaction between fermions with 
opposite spin. We discuss here the case of an infinitely 
long lattice; the finite lattice case with periodic boundary 
conditions (PBC) will be considered elsewhere [l^ . 

There are two similar but different models which leave 
Hamiltonian H invariant under a translation j — > j + 
T (re Z). The first is closely related to the so-called 
Peierls' instability 0, having a r-periodic tunneling rate 
Jj = t + S cos{2nj/T), while its single-particle potential 
Xj vanishes. The second, most commonly used, is that 
of a superlattice, with a constant tunneling Jj — J {> 0) 
and a periodic potential = A cos(27rj7T). We will refer 
to the first of these models as Peierls- Hubbard (PHM), 
and to the latter as ionic Hubbard (IHM). 

The energy bands are obtained by substituting |»I'^) = 

'/'fc,s(j)e"'-'c]<T with 0fc,s(j + t) = 0/c,sO'): into the 
Schrodinger equation "H |^^) = fsl^'a)- For the sim- 
plest case of T = 2, considered from now on, the two 
models have common energy dispersions, namely £s{k) — 
(-1)V(2>/cos(A:))2 + A2, with s = 1,2 labeling the two 
energy bands, and where we have set y/f^ — = J and 
6 = A/2, provided \S\ < \t\. The Bloch functions 4>k,s are, 
on the other hand, different for the two models. Indeed, 
for the PHM <^fc,«(l)/</)fc,,(0) = e''^'--, ^k.s e (0,2^] — 
there is a phase- but not a density-modulation — while 
for the IHM we have 



0fc,s(l) ^ X-Ssjk) ^ 2Jcos(fc) 
(j)k,s{0) 2Jcos(fc) £s{k) + X' 



(2) 



which satisfies |<?!>fe_s(l)| ^ \(j)k,s{0)\, and can even repre- 
sent a particle occupying only one sublattice (odd or even 
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j) for A/J — >■ oo or fc = n/2 (quasi-momenta are defined 
only mod tt). 

The energy dispersions of the PHM and IHM are iden- 
tical and both systems have on-site interaction. Thus, if 
we can analytically solve one of them, the other is solv- 
able as well. In this Communication, we restrict our dis- 
cussion to the IHM since, as the reader can readily check, 
all the formalism discussed below for this model can be 
applied to the PHM. We leave the solution of the PHM 
and further details to a forthcoming paper (l^ . 

c. Wave functions of the two-body problem. Con- 
sider two interacting fermions of different spin compo- 
nents (Ti ^ a2 in the singlet spin state 113| described by 
the IHM with period r = 2. The stationary Schrodinger 
equation for the two-fermion spatial wave function j^*) 
reads, in first quantization, 

-J [*(jl+M,j2)+4'(jl,j2+/i)] (3) 

/J=-l,l 

+ [A((-l)-'"^ + (-Ip) + US,,^,, - E]^{j,,j2) = 0, 

where the wave function is symmetric under the exchange 
of the spatial coordinates, ^'(ji, ^'2) = ^'(j2, ji)- 

In the limit of weak on-site interaction {U = 0), the 
solutions to Eq. ^ are of the form = 

*oO-i,j2) = 65(0fc„.,(ji)0fc„.,(j2)e'(^-^^"^+'=^-'"^)), (4) 

with Os the symmetrization operator [l^ . and eigenen- 
ergies E = Es^^s2-MM = ^siih) + £s2{k2)- For very 
strong interaction \U\/J — > 00 the solutions correspond 
to "fermionized" wave functions [8], 5* = ^i?, 

(jl , J2 ) =6 A (jl )<^fe. (J2)e'('=^^'^ +'-^^'^ ) ) 

xsgn(ji-j2), (5) 

with Oa the antisymmetrization operator (l^ and 
eigenenergies given by £'si,s2;fci,fc2- Here we study the 
more interesting case of finite interaction < \U\/ J < 00, 
for which the analytic solutions were not known previ- 
ously. We discuss below all possible solutions correspond- 
ing to scattering and bound states. 

d. Scattering states. We first discuss the collisional 
states of two fermions. We invoke a periodically modu- 
lated generalization of the two-body Bethe ansatz, writ- 
ten in the most convenient form for our purposes, defined 
for all Jl and j2 as 

^'i3(jl,j2) = *o(ji,i2) + S*fO'i,J2), (6) 

with and 'ip given by Eqs. Q and (O, and B a. c- 
number. The Bethe ansatz ([BJ is an exact eigenfunction, 
i.e. it satisfies Eq. ([3]), if the lattice is homogeneous 
(A = 0). If A is finite, however, for the wave function to 
satisfy at all ji,j2, we need to add an extra wave ^E*!, 
with its corresponding constant C, to the Bethe ansatz. 
The total wave function reads 

«'(ii,j2) = *i3(ji,j2) + C^'iO-i, ^2). (7) 



Clearly, ^'i and must be, asymptotically (|ji — j2\ ^ 
1), solutions to the Schrodinger equation ^ with the 
same energy and total quasi- momentum K — fci + fc2, 
and ^i/^B 7^ const, for aU i? € C. We note that is a 
linearly growing function of — .?2l for vanishing relative 
quasi- momentum, fci = ^2 [l5| ; in such case a generalized 
if-dependent scattering length can be defined as a cx 
[T6| : the divergence of the scattering length denotes a 
bound state entering or exiting the continuum. The cases 
of total quasi-momentum K = 0,7r/2 also have to be 
handled with special care by taking the limits K — > 0, 7r/2 
before using the ansatz ^ in Eq. 

The set of band indices and quasi-momenta 
{si, S2; ^1, ^2} for the Bethe ansatz have the corre- 
sponding energy Esi^s2;ki,k2- Therefore, one has to 
find {s']^, S2; fcj , ^2} so that the following two sets of 
conditions are fulfilled: (li) Es^^s2-MM = £^4,4;fej,fc^ ; 
(lii) k'^ + k2 ~ ki + k2', (2i) if si = S2 = s, either 
s'l 7^ s'2 or s[ — S2 — s and k[ ^ fci,fc2; (2ii) if s\ 7^ S2, 
then s'^ 7^ s'2 and k'^ 7^ ki,k2. Once the primed set is 
found one introduces the wave function ([7]) in Eq. ^ at 
jl = j2 even and odd, thus obtaining a set of two linear 
equations in B and C that can be solved analytically, 
giving in turn the total wave function ^. 

Since we are working with a two-band model, and 
the single-particle bands are mirror symmetric, £s{k) — 
—£s±i{k)^ we have to consider two different possibilities, 
namely intra- and interband scattering. These cases are 
qualitatively different from each other, as we shall see, 
and are therefore considered separately. 

(i) Intraband scattering. We assume that vJ/q^ in- 
terpreted as incident wave, Eq. ([4]), consists of two 
fermions with (real) quasi-momenta fci and fc2, occu- 
pying the same band, si = S2 = s. As discussed in 
the previous paragraph, we have to find a non-trivial 
set {s'l, s'2] k'l, k'2} corresponding to the same energy and 
total quasi-momentum K as the incident state. In the 
majority of the cases, the solution is given by a pair of 
complex quasi-momenta k'l = q = [K + 7r)/2 + iv and 
^2 = 9* (I'ccall K is defined mod tt), with associated 
band indices s'^ — s and s'2 = s ± 1. Such solutions are 
outgoing bound states in the continuum, and represent 
an interband transition [s'^ ^ s'2). Note that interband 
transitions with real quasi-momenta cannot happen since 
we have neither external forces nor energy dissipation 
mechanisms. The only physically acceptable (exponen- 
tially decaying) extra wave with complex quasi-momenta 
is given hy '^1 = '^'^^ , with 

*^r^'0-i-j2) = [^(Ji - j2)</',„.;(ji)</',j,4(j2)e'(^^^"^+^^^-^) 
+ 0{J2 - Ji)0,„4(j2)0,j,4(ji)e'(^^^"=+^=^'^)] 
X e-I^IIJ^-J'il, (8) 

where, in general, qn = in + ^v, and 9 (9) is the step 
function being zero (one) at ji — j2 = 0. In Fig. [T] 
we plot the density j^Pp for a significant case (A/J — 
U/J = 2) in which all the parameters of the Hamiltonian 
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FIG. 1: Density l^*]^ for scattering state with si = S2 = 1, 
fci — 0, k2 = 0.4577. 'I'l corresponds to a bound state in 
the continuum, Eq. (|8]); its contribution is large for small 
interparticle distances (around the diagonal). In this case 
q = -0.27577 + 10.671165, \B\ = 1.915167 and \C\ = 9.719295 
(see text). Values are normalized to the peak. 

are in competition. As observed, the contribution of the 
bound state in the continuum ([8|) is appreciable at small 
interparticle distances, while the wave function behaves, 
for |ji — J2I ^ 1, as the periodically modulated Bethe 
ansatz '^'b- 

This surprising result implies that paired states can be 
created upon two-body collisions on a superlattice even 
without external forces, phonons or dissipation, and are 
due to the combination of the band structure and finite 
interactions. These states should play a non-negligible 
role in the many-electron problem; we have verified that 
the ground state of the system, with PBC, at half- filling 
(4 lattice sites) and quarter-filling (8 sites) is always a 
partially paired state for any non-zero A and U > [12l |. 

There are, on the other hand, cases for which the extra 
wave 5'i has real quasi- momenta k[ and It has now 
the simple form 

^i(ji,j2) = O5(^fc;,4O-i)0fc^,4(j2)e'('='i^'^+'=^-''^)). (9) 

By energy conservation, it is obvious that the associated 
band indices correspond to s[ = s'2 = s. The energies 
and total quasi-momenta for which bound states in the 
continuum or scattered extra waves are needed can be 
inferred from Fig. [51 where the spectrum is plotted for 
X/ J — U/ J — 2. In the figure, the regions of the s — 2 
(s = 1) continuum with energies below (above) the dot- 
ted line (see the discussion on bound states below), cor- 
respond to cases for which the extra wave has real quasi- 
momenta. These regions are already very small - but 
dense - for the value of A under consideration. We note 
that as A/ J gets larger, it becomes less likely that 5'i 
corresponds to a scattered wave. Since, for these cases. 



we can construct two different wave functions (with 
corresponding to ki,k2 or k[,k'2) or any superposition 
thereof having the same total quasi-momentum and en- 
ergy, these states are degenerate (frustrated). 

(ii) Interband scattering. The two incident fermions, 
with quasi-momenta ki and k2 are now in different bands, 
si S2- Energy conservation implies that 4'i has band 
indices s'^ 7^ s'2, and the non-trivial outgoing quasi- 
momenta k'l , k'2 are always real - there is no bound state 
in the continuum - and as a consequence interband scat- 
tering states are always degenerate. The extra wave has 
again the simple form of Eq. and the total wave 

function ^E", Eq. (O, is calculated by substitution in the 
Schrodinger equation at ji — j2, as we have already ex- 
plained. 

e. Bound states. The system under consideration 
also supports bound states - exponentially decaying wave 
functions in the relative coordinate \ji — j2\ with energies 
outside the continuum for each total quasi-momentum K 
- for both attractive and repulsive on-site interaction (if 
U > 0, these are sometimes called antibound states [l3|)- 
Due to the periodicity of the Hamiltonian, we need, as 
for the case of scattering states, two constants so that the 
Schrodinger equation ([3]) can be satisfied for all ji , j2 ■ We 
propose the following ansatz for the wave function 

^{ji,j2)^Ks'Hjuj2) + B^ifHn,j2), (10) 

with Qn = in + ^v, Qn — in + '^v , total quasi- momentum 
= ^1 + ^2 = ^1 + ^2 ,ia], and ^bs given by Eq. ©. 
The two wave functions in Eq. ()10p must have, asymp- 
totically, the same energies, from what we get a relation 
between qn and g^; necessarily, the band indices of the 
functions ^^bs are unequal (except for K — for which 
one of the functions has Si = S2, \v\ < 00, and the 
other has si 7^ S2 with v = 00). The constant B is then 
calculated for ji = j2 = j even (or odd), and qn, (jn are 
varied self-consistently until the Schrodinger equation ^ 
is satisfied for j even and odd, which yields the desired 
bound state energy. Although it is possible to solve this 
problem for general qn and g„, it is simpler to distin- 
guish all possible cases for which Eq. (jlOp is a solution 
to Eq. There are three different types of bound 

states, based on the combinations of qi, qi, as explained 
below. 

Type I: bound states with energies below the Si = S2 = 
1 or above the si = S2 = 2 continuum. These correspond 
to qi = K/2 + iv and qi — [K + 7r)/2 -|- iv, with v ^ v. 
Therefore, the real parts of the quasi-momenta are fixed 
while V and v are calculated self-consistently with B. 

Type II: energies in one of the band gaps and to- 
tal quasi-momenta K G {—ip,ip\, with Q < ip < 7r/2 
depending on the value of U / J and A/J. These have 
qi = {K + 'K)/2-\-iv and qi — [K + -K)/2 + iv, with v ^ v, 
and their calculation is analogous to that of type I. 

Type III: energies in one of the band gaps and total 
quasi-momenta K e (— 7r/2, —ip] U ((/?, 7r/2), with < (/? < 
7r/2 being the same number as for type II if bound states 
of that type exist in their range of quasi-momenta. These 
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FIG. 2: (Color online) Two-body spectrum ior \/ J = U / J — 
2. Shaded regions represent the two-body continua. Black 
dotted lines denote the boundaries between bound states of 
type II and III (see text). Red solid and magenta dashed lines 
are, respectively, energies of states of type II and III. Green 
dashed-dotted lines are energies of bound states of type I. 

have qi — I + iv and qi — £ — iv. One has to calculate, 
for a given starting value of £, the relation between v and 
£ so that the energy is real, and iterate self-consistently 
until the Schrodinger equation is satisfied. 

All bound states for A/J = U/J = 2 are plotted in 
Fig. [21 where we clearly identify three different bound 
state bands. The first one corresponds to a pair bound 
from the lowest scattering band with si = S2 = I, and 
lies in the first band gap. The second band of bound 
states is in the second band gap, and binding energies 
are evidently smaller than those for the first bound band, 
since they are bound from the si ^ S2 band: the density 
modulation prevents high occupancy of two particles at 
the same lattice site if they are in different bands. The 
third band, above the continuum with si = S2 = 2, is 
similar (although not equivalent) to the first bound state 
band. The dotted line in Fig. [5] is the minimum (in 
absolute value) of the energy Ei^2,q,q' with q — {K + 
7t)/2 + iv for each value of the total quasi-momentum K, 
and therefore denotes the boundaries ±ip between bound 



states of type II and III. We note that, for larger values 
of U/ J (not shown), there can be up to four bound states 
at certain values of the total quasi-momentum. 

/. Conclusions. We have found that, on a one- 
dimensional superlattice, the non-trivial underlying pe- 
riodicity has important implications for the two-body 
problem. We have derived exact, analytical solutions 
for the wave functions which show, unambiguously, that 
partial pairing of fermions after their scattering occurs 
over a large range of parameters of the model studied, 
and corresponds to hitherto unexplored bound states in 
the continuum. These states appear only if the incident 
particles occupy the same band. In the case that the 
two fermions occupy different bands, collisions produce a 
phase shift, but also a second outgoing wave correspond- 
ing to quasi-momenta being different from the incident 
ones. This implies that interband scattering is frustrated, 
and the spectrum is degenerate. The implications of our 
results should persist in the solution of the many-electron 
problem at non-zero densities, as can be inferred from an 
analytical two-particle analysis in a finite lattice with pe- 
riodic boundary conditions, which involves a minor gen- 
eralization of our solutions 112]. Our exact results are 
also of relevance to ultracold, low-density gases in opti- 
cal superlattices where the physics is most infiuenced by 
pairwise collisions. Moreover, for bosons, our wave func- 
tions can be used to construct trial functions of p roduct 
(Jastrow) type for the many-body problem [l9| in the 
dilute regime. 
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